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NOTES OR THE PREDICTION OF SHOCK- INDUCED 
BOUNDARY-LAYER SEPARATION 
By Roy H. Lange 


INTRODUCTORY REMARKS 


One of the fundamental problems that appears in the investigation 
of supersonic flow over a surface is that of the phenomena associated 
with the interaction of shock waves and boundary layers. The problem 
of whether a given shock wave will cause bo\mdary -layer separation is 
one which occurs in all cases where a pressure increase is to be obtained 
as a result of the retardation of the flow. Such problems occur, for 
exanqple, in the flow in supersonic diffusers and air inlets and in the 
flow at the rear of airfoils and bodies. Shock-induced boundary-layer 
separation generally results in poor aerodynamic efficiency in the for- 
mer case and in imdesirable airfoil characteristics in the latter case 
and, therefore, this problem is of considerable practical significance. 
The purpose of this paper is to discuss the status of information rela- 
tive to the prediction of shock- induced boundary-layer separation. In 
order to study the fundamental features of the problem, the discussion 
is concerned principally with data obtained on flat plates in two- 
dimensional flow. 


Prandtl has discussed separation of the incompressible boundary 
layer londer the influence of a positive pressure gradient (refs. 1 
and 2). The approximate methods such as those of Von Karman, Pohlhausen, 
and Burl for predicting separation were derived on the assumption that 
the boundary layer has time to adjust itself to a prescribed pressure 
distribution. The Von Karman-Polhausen approximation for a laminar 
boundary layer is 



( 1 ) 


and Burl's approximation for a turbulent boundary layer is 
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( 2 ) 
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boundary- layer thickness 
streamwise pressure gradient 

free-stream dynamic pressure 

en 5 )irical constants 

Reynolds number based on distance 6 

Experience with the use of these approximations has shown that the occur- 
rence of separation depends chiefly upon the pressure gradient dp/dx^ 
and that the turbulent boundary layer can withstand a much greater pres- 
sure Increase before separation than can a laminar boundary layer- When 
the influence of a shock wave on a boundary layer is considered^ it is 
evident that^ if the infinite free-stream pressiire gradient which the 
shock wave represents could extend all the way to the wall^ then separa- 
tion would certainly result; however, as shown in the sketch of figure 1, 
it is known that the pressure difference across the shock is spread out 
in the lower levels of the boundary layer. (See refs. 3 "to 7 ’.) The work 
of Liepmann and Ackeret has shown that the amount of spread of the pres- 
sure rise at the wall depends upon the state of the boundary layer, that 
is, whether the boundary layer is laminar or turbiilent (refs. 3 and 6). 
Thus, the pressure gradient appearing at the wall boundary is fixed by 
the physical properties of the boundary layer and by the strength of the 
shock wave. It seems logical to assume, then, that the occurrence of 
separation in this case depends principally upon the pressure rise P 2 - Pp 

throtigh the shock wave. It was further anticipated that as the pressure 
rise across the shock was decreased there would be one shock strength 
below which no separation of the boundary layer would occur. This con- 
cept was advanced by Beastall and Eggink (ref. 8) and, later, a simpli- 
fied dimensional analysis presented in reference 9 indicated that the 
critical pressixre rise across the shock Z5p/q,p which just causes separa- 
tion of the boundary layer should be proportional to the local skin- 
friction coefficient, c^-. These approximations are extended to the 

case for flat plates in terms of the Reynolds number based on x. Thus, 
for a laminar boundary layer. 


where 

5 

dx 
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and for a turbulent boundary layer 
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It should be en^jhasized that the relationships given in equations (l) to (i^) 
are only approximations. For inconqjressible flow more refined methods have 
been developed (refs. 10 to 15 ) | however, the applications of these methods 
for predicting separation have met with only limited success. A collec- 
tion of the available data for supersonic flow (ref. 9) appealed to bear 
out the predictions shown by equations (5) and (4) at the time they were 
first derivedi however, since that time, more experimental data have come 
to li^t, especially for the turbulent boundary layer, which show that 
the problem must be reexamined. The discussion of these data forms the 
subject of this paper which now follows for both laminar and turbulent 
boundary layers. 


LAMINAR BOUNDARY LAYER 


The available data for shock-induced separation for laminar boundary 
layers on flat plates are given in figure 2 , where the critical pressure 
rise Ap across the shock divided by the free-stream dynamic pressure qj^ 
is plotted against Reynolds number on logarithmic scales. The Reynolds 
number is based on the distance from the leading edge of the plate to 
the point of intersection of the shock wave 'and the boundary layer. The 
sources of these data are given at the top of the figure. (See refs. 5j 
6, 8, II 4 -, 15> and 16 .) The data at Mach numbers of 1.4, 1.44-, 2.00, and 
2.05 (refs. 3> 16 , and 15 ) were obtained from tests in which shock waves 
of varying strength were made to in^iinge upon the boundary layer on a 
flat plate. The data at Mach numbers of 1.93 and 2.48 (refs. l4 and 8 , 
respectively) were obtained in the separated region ahead of a forward- 
facing step. It can be seen that the available data are rather limited 
in scope arid, therefore, are not conclusive j however, there are some 
trends in the data which should be mentioned. For example, at free- 
stream Mach nimabers Mp of 1.93? 2.00, 2 .O 5 , and 2.48 the Reynolds nirni- 
ber effect on the critical pressure coefficient appears to follow the 
inverse square root of the Reynolds number as denoted by the dashed lines 
in figure 2. Except for the data at Mach numbers of 1.40 and 1.44, the 
critical pressure coefficient also decreases with Increasing Mach number. 
These trends of Reynolds number and Mach niimber agree with the predictions 
of equation (3)j however, the magnitude of the Mach number effect shown, 
especially between Mach nimbers of 1 and 2, is much greater than that 
which would be predicted by reference 9> Recent data obtained at the 
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Ames Laboratory in the separated region aiiead of a forward-facing step 
show an increase in Ap/q^^ with increase in Reynolds mimber which is 
opposite to that obtained at Mach numbers of 1«95 and 2.,hQ. The forward- 
facing-step data shown in figure 2 should be ignored, therefore, until 
more systematic data are available. Stewartson (ref. I7) has made a 
detailed analysis of the interaction process which leads to the inference 
that the dimensionless pressure rise required to produce laminar separa- 
tion would be proportional to Also shown in figure 2 is a curve 

which traces the criterion of separation advanced by Pabst (ref. 18) in a 
recent Argentine paper| however, this criterion cannot account for the 
Mach number effect and does not correlate with any of the experimental 
data shown. 


TURBULEMT BOUNDARY LAYER 


Investigations of shock—bovindary-layer interaction for the t-urbu- 
lent boundary layer have shown that a given shock wave may or may not 
separate the boimdary layer.- Data are now available from a number of 
sources in which turbulent boundary -layer separation has been investi- 
gated by three methods; (l) the forward-faclng-step technique, (2) the 
wedge technique, and (3) the incident -shock technique. 

In order to remove all doubt as to whether the turbulent boundary 
layer has been separated, several investigators have forced separation 
by means of a forward-facing step mounted on a flat plate (see refs. 8, 

9, 19 and 20). Typical data for this type of configuration are given 
in figure 3 which shows the pressure distribution along the surface and 
(to the same scale) a sketch of the flow field in the interaction region 
as determined from shadowgraphs. These data were obtained in a blowdown 
jet of the Langley gas dynamics laboratory at a Mach number of 5-05- The 
flow diagram at the top of the figure shows that a wedge-shaped separa- 
tion region is formed ahead of the step and is bounded on its upstream 
edge by the shock wave. The direction of the circulatory flow within 
the separated region is shown by the arrows. 

The pressure coefficients on the plate first reach a maximum value, 
noted herein as the first peak, at a point about halfway between the 
location of the shock wave and the location of the step. This distance 
is roughly the equivalent of 8 boundary-layer thicknesses or 133 momenttim 
thicknesses, on the assixngition of a l/7-power velocity distribution in 
the boundary layer just ahead of the shock. The pressures then dip 
slightly behind the first peak and siibsequently rise sharply, showing 
the large influence of the circiilatory flow. Also pertinent to the dis- 
cussion of the flow in the separated region are the pressure coefficients 
measured along the front vertical face of the step given in figure k. 
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The three isolated points at a Reynolds mamher of 4 X 10^ were obtained 
at Mj_ = 1,86 (ref. 2l) ^ and the data for Reynolds numbers ranging from 

12 X 10^ to 32 X 10^ were obtained at = 3*03* The pressure orifices 
were located at the base of the step and at two other vertical locations 
above the surface of the plate as denoted by z/h» The data at Mp = 3*03 
show no significant Reynolds mmiber effect on the pressure coefficients. 

The results show that there is one stagnation point at the foot of the 
step and one near the top of the step^ and calculations based on the data 
at Mp = 3*03 and utilizing the incompressible Bernoulli equation show 
that the velocity downward along the vertical face is about one-fourth 
the free- stream velocity | whereas the velocity along the plate in a 
direction opposite to the main flow is about one-third the free-stream 
velocity. Thus^ the separated region cannot be treated as a dead-air 
space as is commonly assumed. The results at both Mach numbers also 
show that a considerable error would result if the pressures on the 
front face of the step were assumed to be the same as that obtained on 
the plate surface ahead of the step in the separated region. The first 
peak pressure coefficients obtained ahead of the step are shown by the 
dashed lines at both Mach numbers for comparative purposes in this case. 

It is clear then^ from the results given in figures 3 sjid that the 
first peak pressure coefficient is obtained as a result of the mutual 
effects of the shock on the boundary layer and of the circulatory flow 
in the separated region and should not be interpreted as the value of the 
pressure rise across the minimum strength of shock wave which just causes 
separation of the boundary layer. These results have been obtained for 
cases where the step height is about 3 times the local boundary-layer 
thickness and may be changed somewhat for cases where the step height 
is very large compared with 6. 

A summary of the available data obtained from the use of the step 
technique for forcing boundary-layer separation is given in figure 5 
which shows taken at the first peak plotted against Reynolds num- 

ber on a logarithmic scale* The Reynolds number is based on the distance 
from the leading edge of the plate to the point of intersection of the shock 
wave with the boundary layer. All the data were obtained from pressure 
distributions (see refs. 8^ 1^;, 20^ 22 j, and 23 )^ and the sources are 
given at the top of the figure. The Mach number range of the data is 
from 1.55 shown by the long string of points at the top of the data to 
3»65 shown by the lowest data points « The pressure- distribution data 
at M]_ = 3*03 given by the circles are new data which have not been^ 

published. The data given in reference 9 (NACA TN 2770 ) for Mj_ = 3*03 
represented by the dashed line which varies as Rx“^/5 were obtained 

by measuring shock angles close to the point of intersection of the 
shock wave and the boundary layer ^ where ^ as shown previously^ the pres- 
sures on the plate are changing rapidly ^ therefore this method for 
obtaining pressure coefficients is too crude and the data should be 
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ignored. It is apparent from the mass of data that, except for the data 
at Mach numbers of 1.86 and 2.48, the Reynolds number effect on the value 
of (Ap/qj^) 2 ^g^ peak very slight. On the basis that there is no 
Reynolds niimber effect, figure 6 has been prepared to show the decrease 
in (Ap/q.p)p 3 t peak increase in free-stream Mach number for Mach 

numbers between 1-55 and 3 . 65 . All the data from the previous figure 
have been included in this plot, and the vertical lines connecting some 
of the symbols show the extent of the Reynolds number effect obtained. 
Included on this plot is the empirical relationship derived by Beastall 
and Eggink from a curve which best fit their data for both forward- 
facing steps and backward-facing steps (refs. 8 and 24). This approxi- 
mation is independent of both Reynolds nunber and Mach nunber and, there- 
fore, does not correlate well with the aArallable experimental data for 
forward-facing steps. 


The second technique for producing tvirbulent boundary-layer separa- 
tion is the use of wedges of different angles moionted on flat plates, 
and a limited amount of data is available, (See refs. 20 and 25.) This 
configuration is analogous to the deflection of a flap or a control sur- 
face. Typical data obtained at a Mach nunber of 3-03 are given in fig- 
Tore 7 which shows the pressure distribution along the plate and on the 
wedge and above it a sketch of the flow phenomena as determined by shadow- 
graphs. A double scale is given along the abscissa of the pressure dis- 
tribution - one which gives x in Inches measured from the leading edge 
of the wedge and one which gives a measure of the botmdary-layer thick- 
ness, x/ 6 . As shown in the flow picture, the separation in the corner 
produced by this particular wedge angle results in a weak shock wave, 
which projects ahead of the main shock, and an inflection point is obtainec 
in the presstire distribution on the surface. Downstream of this point 
the pressure coefficient continues to rise and levels off at a value 
somewhat less than that calculated from oblique-shock theory for this 
wedge angle in the absence of a boundary layer. In general, the limited 
available data at a given Mach number show that, for wedge angles greater 
than a certain value, the pressure distribution has an inflection point 
similar to that shown in figirre 7 j moreover, the value of /^p/qp meas- 
ured at the inflection point remains almost constant with f\irther increases 
in wedge angle. The data at Mp = 3»03 also show that the value of Z^p/qj 

obtained at the inflection point is essentially constant for Reynolds num- 
bers ranging from 12 X 10^ to 32 X 10^. Results are available from tests 
utilizing the third technique in which shock waves of varying strength 
are made to inpinge upon the boundary layer on a flat plate. (See 
refs. 16 and 26 .) In these tests inflection points are obtained in the 
pressure distributions along the plate swface somewhat similar to those 
in the wedge tests, and these inflection points are also associated with 
local separation of the turbulent boundary layer. The tests of Gadd and 
Holder at a Mach number of 2 show no significant effect of Re 3 molds num- 
ber on the value of /p/qp obtained at the Inflection point for Reynolds 
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numbers ranging from atout 0.8 x 10^ to 10 x 10^. In figure 8 Ap/q^ 

is plotted against Mach number, where the inflection-point pressure 
coefficients obtained in the wedge tests are given by the open symbols 
and the inflection-point pressure coefficients obtained by the incident- 
shock technique are given by the solid syribols. Also shown on this fig- 
ure is the curve representing the data obtained by the forward-facing- 
step technique. The data given on this figure, therefore, constitute 
all information available at present on turbulent boundary-layer separa- 
tion. The spread in /ip/q^ obtained at = 1.80 in the wedge tests 

represents a Reynolds number effect, although, as mentioned previously, 
no such Reynolds number effect was obtained at = 5>03* The spread 

in /^P/qq at = 2 in the incident-shock tests represents the maxi- 
mum scatter in the data. Althoiigh the available data are rather limited 
in scope, the results show that the inflection-point pressure coefficients 
obtained from both techniques generally have the same range of values 
with increasing kfeich number and that on the average these values are 
about 20 percent lower than those obtained using the step technique. 

The application of these data for predicting separation shotild, there- 
fore, be limited to these particular configurations, at least for the 
present. For example, the data from the incident-shock technique repre- 
sent conditions of local separation of the flow and, because the experi- 
ments are performed on flat plates, the flow reattaches downstream of the 
separation point. This reattachment may be changed somewhat for condi- 
tions where a back pressure exists - for example, for conditions near the 
trailing edge of an airfoil. Also, flight data for a wing in transonic 
flow indicate that the /!p/qq for separation is predicted more accurately 

by the step data if extrapolated to the lower supersonic Mach numbers 
obtained in the flight tests (ref. 27). These data are usefiiL, then, in 
providing a first approximation to the pressure coefficient for which 
separation is likely to be encountered. 


CONCLUDING REMARKS 


In conclusion, the present status of information relative to the 
prediction of shock-induced boundary-layer separation indicates that, 
althou^ no universal value of pressure-rise coefficient which causes 
incipient separation of the bomdary layer has been found, there is a 
fairly narrow band of pressure coefficients from which predictions of 
turbulent separation can be made with an accuracy probably sufficient 
for engineering purposes. On the basis of these results the following 
tentative conclusions are given; 

1. The data obtained with forward -facing steps, wedges, and inci- 
dent shock waves indicate that there is a dependency of the pressure 


CONFIDENTIAL 



8 


CONFIDENTIAL 


NACA RM L55Il6a 


coefficient for separation on Reynolds number for the laminar boundary 
layer but little, if any, dependency on Reynolds nimiber for the t\irhu~ 
lent boundary layer. There is a dependency of this pressure coefficient 
on Mach number for both laminar and turbulent boundary layers. 

2. For the particular case of the spoiler, the available data 
obtained by the forward -facing-step technique permit calculations of 
the loading on the sixrface ahead of the spoiler, the pressure on the 
front face of the spoiler, and the separation point ahead of the spoiler 
for a Mach number range of from 1.55 'to 5*65 for the turbulent boundary 
layer. 

5 . For application to supersonic diffusers or scoop Inlets, the 
available data from Incident-shock-wave tests provide a first approxi- 
mation to the minimum strength of shock which will separate the turbu- 
lent boundary layer for Mach numbers between 2 and 3- 

h. From the data available from the wedge tests, a first approxi- 
mation to the pressure coefficient for which separation becomes apprecia- 
ble as a result of flap deflection can be made for a surface with a tiurbu- 
lent boundary layer for Mach numbers between 1.75 and 5»03. 

5. Caution should be exercised in atten^jting to predict the separa- 
tion or loading on configurations which differ considerably from those 
for which experimental data are available. For example, fair success 
has been obtained in predicting base pressirre coefficients by the use 
of the forward-facing-step data, but reasons for this success are not 
at present fully understood. 
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SHOCK WAVE ENTERING A BOUNDARY LAYER 
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LAMINAR BOUNDARY LAYERS 
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SEPARATION OF TURBULENT BOUNDARY LAYER BY STEP 
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FIRST PEAK PRESSURE COEFFICIENTS AHEAD OF STEPS FOR 
TURBULENT BOUNDARY LAYERS 
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Figure 5 
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ABSTRACT 


The present status of available information relative to the predic- 
tion of shock- induced boundary- layer separation is discussed. Experi- 
mental results showing the effects of Reynolds number and Mach number 
on the separation of both laminar and turbulent boundary layer are given 
and con^jared with available methods for predicting separation. The flow 
phenomena associated with separation caused by forward-facing steps, 
wedges, and Incident shock waves are discussed. Applications of the 
flat-plate data to problems of separation on spoilers, diffusers, and 
scoop Inlets are indicated for turbulent boundary layers. 
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